Quantum impurity spin in Majorana edge fermions 
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We show that Majorana edge modes of two-dimensional spin-triplet topological superconductors 
(superfluids) have Ising-like spin density whose direction is determined by the d-vector characterizing 
the spin-triplet pairing symmetry. Exchange coupling between an impurity spin [S = ^) and 
Majorana edge modes is thus Ising-type. Under external magnetic field perpendicular to the Ising 
axis, the system can be mapped to a two-level system with Ohmic dissipation, which is equivalent 
to the anisotropic Kondo model. The magnetic response of the impurity spin can serve as a local 
experimental probe for the order parameter. 

PACS numbers: 



Majorana fermions are fermionic particles that are 
their own antiparticle. Originally proposed long ago 
to describe neutrinos in high energy physics Majo- 
rana fermions have recently been a subject of intensive 
studies in condensed matter physicsi^i^ Their mixed na- 
ture of being particle and antiparticle implies that Ma- 
jorana fermions may emerge as elementary excitations 
in superconductors and superfluids where the number 
of particles is not well-defined. Indeed they are theo- 
retically predicted to appear as gapless boundary exci- 
tations of topological superconductors/superfluids with 
spin triplet pairing)^ Candidates of such topological 
materials include superfluid phases of '^Hcj^ the su- 
perconducting states of Sr2Ru04^ and possibly some 
non-centrosymmetric superconductors.^ However, being 
"real-part" of ordinary (complex) fermions, Majorana 
fermions are charge neutral and only weakly interacting 
with other particles; they are hard to detect. To probe 
and control Majorana fermions is therefore a great chal- 
lenge. 



In this paper we study a quantum impurity spin cou- 
pled with Majorana edge modes with spin degree of free- 
dom. We first argue that a distinct signature of Majorana 
edge fermions of a generic two-dimensional (2D) spin- 
triplet topological superconductor (superfluid) is Ising 
character of their spin densityj^iiS To probe this Ising 
spin, we consider a spin-i magnetic impurity coupled 
to the Majorana edge modes. We show that the impu- 
rity spin has strongly anisotropic and singular magnetic 
response which is due to quantum dissipation from the 
Majorana edge modes. We propose that electron spin 
resonance (ESR) can serve as a novel local probe for the 
Majorana fermions in spin-triplet topological supercon- 
ductors and superfluids. 



Let us take a closer look at the Majorana edge modes 
of 2D spin-triplet topological superconductors. We first 
consider the BdG Hamiltonian of a prototypical topolog- 
ical superconductor, a spinless chiral p-wave supercon- 



ductor with px ± ipy symmetry 
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V 2ikF 



(1) 



Here r — {x,y), d± — ± idy, d = 9+9-, kp the 
Fermi wave number, ^ — h kp/2m, and is a U(l) 
phase. In the bulk superconductor with a spatially uni- 
form pair potential A(r) = A, quasiparticle spectrum is 
fully gapped. When the superconductor has a boundary, 
the BdG Hamiltonian has a gapless edge mode with lin- 
ear energy dispersion ^ vk — /S.k/kp {\k\ < kp). For 
a straight boundary defined hy X = x cos 4>-\-y sin = 0, 
we solve the BdG Hamiltonian on the half plane X < Q 
[A(r) = A0(— X)], assuming fixed boundary conditions 
aX X — Q. We then obtain the edge mode's wave function 



Uk{r) 
Vk{r) 



Wk{X) 



„i(2i5±20+7r)/4 
-i(2)S±20+7r)/4 



(2) 



Here Y is the coordinate along the edge, Y = —xsm(j)-\- 
ycoscj), and Wk{X) is the normalized real- valued wave 
function localized at the edge. The particle-hole sym- 
metry, VHj-V'' ~ —Hj- with V — aiK, guarantees that 
{vl,u^) = {u^k,V-k) is also an eigenmode with energy 
—vk, where K is complex conjugation and is the jth 
component of the Pauli matrices. The mode expansion 
of the Nambu field (-0,^'^)* for l-^'l < A is then given by 



0t(r) 



dk 



ul{r) 



(3) 



which leads to the following condition of Majorana type 



V'(r) 



(4) 



There are two distinct types of 2D spin-triplet topo- 
logical superconductors (a) chiral type^ and (b) he- 
lical type.— '^^1^^ They can be obtained, for example, by 
combining two copies of the spinless chiral p-wave su- 
perconductors of same or opposite chiralities; their BdG 
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Hamiltonians are given by Ti.^ Ti.^ and Ti.^ © T-Lj; , re- 
spectively. The two types of 2D spin-triplet topological 
superconductors are characterized by the order parame- 
ter Afc = idk ■ (T(J2 with the d-vectorjii 



{xs0 + ycg)ik.^ 
{xsg + ycg)kx -( 



iky), (chiral), , . 

{xcg — ysg)ky, (helical). ^ ' 



Here {cg,sg) = (cos 61, sin 0), 9 — ^(6*^ — 6*^), and x and 
y are unit vectors in the spin space. Notice that the 
direction of d^ is independent of k in the chiral case, 
while it constitutes a coplanar spin texture in the fe-space 
in the helical case. 

The two types of 2D spin-triplet topological super- 
conductors support gapless Majorana edge modes, each 
spin component of which obeys Eq. ([4]): (a) ipair) = 
"^■01. (r) (fT =1 and I) for the chiral case, and (b) 



le 



^^{r) = ie*^t+#^t(^) ^ ie'^i-^'i'il)'^^{r) for the 

helical case. These Majorana conditions lead to the op- 
erator identities for the edge mode's spin density. 



,-2i(e+0)- 



(chiral). 



s_(r) (helical). 



(6a) 
(6b) 



These equations imply that the spin density is always 
Ising-like: {sx,Sy) cx {sg,cg) for the chiral case and 
{sx,Sy) oc (s0+0, C0_|-0) for the helical case. Comparing 
Eq. dS]) with Eq. (l6b)) . we can see that this Ising direc- 
tion is dictated by the d-vector in the bulk as. 



s(r) cx 



dfc (chiral), 
dj^ (helical). 



(7) 



where X = (cos 0, sin 0) is a vector normal to the bound- 
ary. In the helical case, s(r) depends on the direction of 
the boundary, as the rotation in the real space is tran- 
scribed into that in the spin space. Otherwise, it does 
not depend on the shape of the boundary. 

The Ising-like spin density ([7|) is a hallmark of the 
Majorana edge modes. Its strong anisotropy reflects the 
spin-triplet pairing symmetry in the bulk topological su- 
perconducting order. In the following we will study quan- 
tum dynamics of a magnetic impurity (probe spin) inter- 
acting with this Ising spin density. 

The Hamiltonian for the coupling between the Majo- 
rana Ising spin density s{r) and a spin-i probe at r = 0, 
S = [Sx, Sy, S^), is given by 



(8) 



Here the probe spin's direction is taken to be paral- 
lel to s{0), and ipcr{Y) is a one-dimensional (ID) Majo- 
rana field satisfying = ■0^ and {'tpa{y),4'l,,{Y')} — 
Sa,<j'HY — Y')- The ID Majorana field is obtained from 
ipaif) by appropriate U(l) gauge transformation and 
dropping unimportant X-dependence [i.e., Wk{X) 1]. 



The Kondo coupling in Eq. 
Anderson impurity model. 



can be obtained from the 



4^,(0) + h.c.]} + 



with rid, (J = d\dcr, by the standard procedure. This yields 

= -i di[ay\c,pdp and J = 2eUl[{U + ed)\ed\] (> 0). 
The kinetic energy of the Majorana edge modes reads 



"Hkin = iv 



dY 



(9) 



where the +/— signs in the integrand are for the chi- 
ral/helical superconductors, respectively. The ground 
state of 'Hkm + "Hex is doubly degenerate, Sz = ±i. 

We will make full use of the knowledge from ear- 
lier studies on dissipative two-state systems and Kondo 
effect,—^— to show that the impurity spin, when sub- 
jected to external magnetic fields, exhibits anisotropic 
dissipative quantum dynamics due to the "background" 
Majorana edge modes. 

Consider first that a magnetic field is applied perpen- 
dicular to the Majorana Ising spin, say along the x- 
direction; the Hamiltonian reads 



n = n 



kin 



He 



hSx 



(10) 



Interestingly, the system undergoes a quantum phase 
transition, depending the exchange coupling J (Fig. 1). 
To see this, let us map the Hamiltonian onto the Ohmic 
dissipative two-state system,— 

{dy'l>fdy+^S4dy'P)\y = + hSx, (H) 



/2tt 



where the bosonic field $(a;) obeys [$(x), 9y$(y)] = 
i5{x — y). We have combined the two species of chi- 
ral Majorana fields, {ip^,ipi}, into a single spinless chi- 
ral fcrmion, ^'(y) = [ipf{y) + iipi{±y)]/V2, where the 
+/— sign refers to the chiral/helical superconductors, 
respectively. The Ising spin density at r = is then 
reduced to the fermion number density 4'^(0)\E'(0). The 
complex fermionic field is readily bosonized in terms of 
the phase operator, ^'(y) oc exp[i\/27r$ (?/)], which trans- 
forms Eq. (fTO|) into Eq. The Ising exchange cou- 
pling J and the transverse field h control the coupling 
to the Ohmic bath and tunneling between the two states 
{Sz — 5, — ^), respectively. The dissipative two-state 
system is related to the anisotropic Kondo modelj ^^i^^ 

Hk - -2tv I dy [f\{y)dy'f^{y) + vl^|(y)9,vl;^(y) 

+ Ji_ [5+*|(0)«'t(0) + 5_*|(0)«';(0) 

+ JzSz [*|(0)*t(0) - «'l(0)«'i(0) 



(12) 



where ^'-1-4 are complex spinful fermionic fields and 5-1- 
Sx iiSy. The coupling constants are related by 



^TTV 2tTV \ 2'KV 



(13) 
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Toulouse limit 
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FIG. 1: (color online) Renormalization group flows. 



where r is inverse of an ultraviolet cutofF (band width) of 
the massless Majorana edge modes, i.e., r^^ = A. The 
transverse (J_l) and longitudinal (Jj) Kondo exchange 
couplings are related to the transverse field h and the 
Ising exchange coupling J, respectively. With the help 
of these mappings, the renormalization group (RG) equa- 



tions for rj = hr and e 
tainedfi^ 



(J/4. 



can be readily ob- 



de 
dlnr 



= -W, 



These scaling equations lead to the RG flow diagram in 
Fig. m which has a fixed point at (J, h) = (0,oo) and a 
line of fixed points at /i = for | J| > 4\/27rt; = Jc (e > 2). 
When I J| < Jc, the transverse field his a. relevant pertur- 
bation, and the system flows toward the high-field fixed 
point. In this case the ground state is a quantum me- 
chanical superposition of Sz — +\ and — ^ states, corre- 
sponding to the antiferroniagnetic Kondo singlet. When 
I J| > Jc, on the other hand, the orthogonality catastro- 
phe of Majorana edge excitations suppresses the tunnel- 
ing, thereby making a small transverse field h irrelevant. 
For small fields h < he, where 



(15) 



the RG flows end at the zero-fleld fixed line; this corre- 
sponds to the ferromagnetic Kondo regime, Jz < — | Jj.|. 
For large enough fields h > he, the system is eventually 
renormalized to the strong-coupling regime (see Fig. [ij. 
The crossover from weak- to strong-field regime occurs 
around the Kondo temperature Tk given asr^ii^ 



_ j Aexp(-6/\/^ 



'), |J| > Jc n /i > he 

\J\ < Jc, 



(16) 



where h — tt A / V8/ic- 

Let us assume that the magnetic field has a component 
in the direction of the Majorana Ising spin. This adds 
a weak Zeeman term T-Lz = —hzSz to the Hamiltonian. 



From the mapping to the anisotropic Kondo modelj^^— 
we see that the response of the probe spin S to a weak 
longitudinal field hzZ should be singular at the critical 
field h = he for | J| > Jc. In the ferromagnetic Kondo 
regime where h < he and |J| > Jc, the ground state 
is doubly degenerate, and the magnetization Mz = {Sz) 
shows a jump across /i^ = at zero temperature As h 
is increased, the magnetization jump is reduced but re- 
mains finite at h = he<^ At finite temperature T, the 
static longitudinal susceptibility Xzz = dMz/dhz\h^=o 
shows a Curie behavior, Xzz oc 1/T, as long as ft. < he- 
Above the critical field he, the ground state is a singlet. 
In this case, the longitudinal susceptibility saturates be- 
low the Kondo temperature, Xzz ~ ^/Tk at T < Tk- 
The saturated value diverges at h = he- 
Ac longitudinal field, hzZe~^'^^ , induces transitions 
between the levels split by the tunneling term hSx- 
To see this, let us consider the dynamical susceptibil- 
ity x..(w) = I /o°°dte'(-+«)*([&(t),^,(0)]), which is a 
Fourier transform of the time-evolution of the magne- 
tization, P{t) = {Sz{t)), studied in the context of the 
dissipative two-state system.— The transition spec- 
trum is obtained from the imaginary part of Xzz{^) 
and has a universal scaling form which depends on e, 
Szz{^) ■— ImXzz('^)/w = fe{uj/TK)- Its qualitative fea- 
ture is well understood (for small h/A),- 



16.17.19 



as we 



(14) briefly summarize below. 



In the weakly dissipative regime e < |, Szzi^^) has a 
peak at w « Tk, which signifies coherent transitions be- 
tween the bonding and anti-bonding states split by the 
renormalized transverse field (Fig. [^b)- Stronger dissi- 
pation destroys the coherent tunneling; when e > |, the 
coherence is lost and 6*22 (w) shows only a diffusive peak 
at w = (Fig. [2Jl). The half- value width of this peak 
is on the order of Tk-^^ As the dissipation is further in- 
creased toward the ferromagnetic Kondo region, the dif- 
fusive peak gets sharper, as Tk — > at e = 2. 

In the ferromagnetic Kondo region, the spectral func- 
tion at T = has a (5- function peak at w = OJ^^^ 
At finite temperature, or in the presence of a finite 
dc "bias" field hz, the width of the diffusive peak at 
a; is finite, as T and hz play a role of low-energy 
cutoff.iS'i"^ The half-value width ujQ can be obtained by 
2nd order perturbation in the transverse field h, yield- 
ing ujo (X (/iVA) (T/A)^-i or {h^ / A){hz/ Ay-\ A 
diffusive peak with the width of the same T and /in- 
dependences also appears in the antiferromagnetic Kondo 
regime, when T > Tk or hz > Tk- 

Next, we discuss response of the probe spin S to 
the transverse field h (while we set hz = 0), such as 
the magnetization — (Sx) and the susceptibility 
Xxx = dMx/dh- This response is unique to our system 
and has not been considered in the related two models 
(ITT]) and ((HI), as the Zeeman term hSx corresponds to 
the tunneling and the Kondo exchange interaction, re- 
spectively. 

At the critical field h = he separating the ferromag- 
netic and antiferromagnetic Kondo phases, the transverse 
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edge singularity, 



20 



FIG. 2: (color online) Transition (absorption) spectra, (a) 
S± (iu) = Im Xxx (i^) /i^ under longitudinal dc field and weak 
dissipation e < 1. (b) 511(0;) = ImXzzi^)/^^ under transverse 
dc field h and weak dissipation e < |. (c) S±{u}) for e > 1, 
and (d) Sy (tj) for e > | {Tk for e > 2). 



magnetization shows only weak (essential) singular- 
ity, Mx{h) - M^{hc) oc exp(-6/^|/i| - he). This sin- 
gularity is like the specific heat anomaly at Kosterlitz- 
Thouless transition and is too weak to be observed. Thus 
Mx increases monotonically as a function of h without 
any hint of anomaly at h — he- 

In the limit of small transverse field h ^ A, the trans- 
verse magnetization increases linearly with h, when 
dissipation is strong, e > 1. The susceptibility Xxx is di- 
verging as (e — 1)^^ for e 1, and at the point e = 1, the 
linear h dependence acquires a logarithmic correction, 
Mx = ln(/i/A) -I- • • • . This result can be ob- 

tained from the exact ground state energy in the Toulouse 
limit In the weakly dissipative case < e < 1, the h 
dependence becomes sublinear, oc (/i/A)'^^*-^"'-'. This 
indicates that the linear susceptibility diverges in the low- 
T limit. However, this singularity is weaker than that of 
the Curie behavior and is given by 



sm 



ttT\ 



l(xT-l+^(17) 



In the presence of dc longitudinal field hzZ, a small ac 
transverse field /lie^*"* induces transitions between the 
Sz — ±5 states. The transition spectrum has a divergent 



Imxx£c(ii^)| 



T=h=0 



f(^A^^737i:Re(--M,(i8) 



for positive uj, where c and c' are some positive constants. 
Equation (ITS)) is also valid for e > 1. 

Based on the results discussed so far, we propose that 
ESR measurements of the impurity spin may be used to 
identify the direction of the Majorana Ising spin. A dc 
field smaller than the (lower) critical field i?c(i) will cause 
precession of only those probe spins located around the 
boundary (edge). An additional ac field (perpendicular 
to the dc field) will then induce resonant transitions. Fig- 
ure O shows qualitative picture of the absorption spectra 
for two complementary experimental geometries, i.e., the 
dc field applied parallel (a,c) and perpendicular (b,d) to 
the Majorana Ising spin. The upper two panels (a,b) as- 
sume weak dissipation. The spectra show the strongest 
edge singularity at the Larmor frequency hz, when the 
dc field is parallel to the Majorana Ising spin (a). When 
the spin precession is driven by the transverse field h 
(b), the singularity is replaced by a resonance peak at 
Lj = Tk = /i(/i/A)'/(2-') < h. The clear difference be- 
tween the spectra measured in the two experimental ge- 
ometries persists at stronger dissipation (compare c and 
d) . We thus expect that the Ising direction can be identi- 
fied in principle by measuring the absorption spectra for 
various directions of the dc field. Incidentally, anisotropy 
in the static susceptibility may also serve as a good indi- 
cator; it shows Curie-law along the Ising direction (xzz), 
while it remains constant for e < 1 or diverges as T'^~'^ 
for e > 1 in the perpendicular direction {xxx = Xyy)- 

In summary, we proposed to probe the massless Majo- 
rana edge modes of spin-triplet topological superconduc- 
tors and superfluids, by introducing a quantum impurity 
spin to their boundary. Anisotropic magnetic responses 
of the probe spin directly reflect the Majorana nature of 
the edge excitations. 
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